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Gravitoelectromagnetism is briefly reviewed and some recent developments in this 
topic are discussed. The stress-energy content of the gravitoelectromagnetic field 
is described from different standpoints. In particular, the gravitational Poynting 
flux is analyzed and it is shown that there exists a steady flow of gravitational 
energy circulating around a rotating mass. 



1 Introduction 

Gravitoelectromagnetism (GEM) is based upon the close formal analogy be- 
tween Newton's law of gravitation and Coulomb's law of electricity. The New- 
tonian theory of gravitation may thus be interpreted in terms of a gravito- 
electric field. Any field theory that would bring Newtonian gravitation and 
Lorentz invariance together in a consistent framework would necessarily con- 
tain a gravitomagnetic field as well. In general relativity, the non-Newtonian 
gravitomagentic field is due to mass current and has interesting physical prop- 
erties that are now becoming amenable to experimental observation. 

Developments in electrodynamics in the second half of the nineteenth cen- 
tury led Holzmiiller [1] and Tisserand [2] to postulate a gravitomagnetic com- 
ponent for the gravitational influence of the Sun on the motion of planets. The 
magnitude of this ad hoc component could be adjusted so as to account for 
the excess perihelion motion of Mercury. However, Einstein's general relativ- 
ity successfully accounted for the perihelion precession of Mercury by means 
of a post-Newtonian correction to the gravitoelectric influence of the Sun. The 
general relativistic effect of the rotation of the Sun on planetary orbits was first 
calculated by de Sitter [3] and later more generally by Thirring and Lense [4] . 
Evidence for the Lense-Thirring precession of laser-ranged satellites LAGEOS 
and LAGEOS II around the rotating Earth has been described by Ciufolini [5]. 
The rotation of the Earth generates a dipolar gravitomagnetic field; a direct 
measurement of this field via the precession of gyroscopes in orbit about the 
Earth is one of the main objectives of NASA's GP-B that will be launched in 
the next few years [6]. 

The Lense-Thirring orbital precession and the gravitomagnetic gyroscope 
precession point to a certain temporal structure around a rotating mass. This 
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is best expressed via the gravitomagnetic clock effect [7]. In it simplest form, 
the clock effect involves the difference in the periods of two clocks moving in 
opposite directions on the same orbit. Let r+(r_) be the proper period of a 
standard clock that makes a complete revolution around a rotating mass on a 
prograde (retrograde) orbit; then, for a circular orbit of radius r ^> 2GM/c 2 
in the equatorial plane, t + — t_ w AttJ/Mc 2 . The gravitomagnetic clock 
effect is in the lowest order independent of Newton's gravitational constant 
and the radius of the orbit. The key requirement for the existence of this effect 
is azimuthal closure [8]. For a circular equatorial orbit around the Earth, 
t + — t_ « 10~ 7 sec; however, the detection of this effect would involve a 
difficult experiment [9]. Various complications have been discussed in recent 
studies [10]. An important feature of the gravitomagnetic clock effect is that 
the prograde period is longer that the retrograde period; this is contrary to 
what would be expected on the basis of the heuristic notion of "dragging" 
associated with rotating masses [8]. 

There are different approaches to GEM within the framework of general 
relativity [11]. For the sake of simplicity and convenience, we adopt a conven- 
tion that results in as close a connection with the formulas of the standard 
electrodynamics as possible [12]. In this convention, a test particle of iner- 
tial mass m has gravitoelectric charge q^ = —m and gravitomagnetic charge 
qg = —2m. If the source is a rotating body of incrtial mass M, the corre- 
sponding gravitational charges are positive, i.e. Qe = M and Qb = 2M, in 
order to ensure that gravity is attractive. Thus, we always have qs/qE = 2; 
this can be traced back to the spin-2 character of the gravitational field. Hence 
for a spin-1 field qB/qE = 1, as in Maxwell's theory. 

2 GEM Field Equations 

Let us begin with a linear perturbation of Minkowski spacetime given by g^ v = 
V/j-v + hfxu, where = (ci,x). It is useful to define the quantity h^ v — 
hfxu — \ fliLvh, where h =tr(/i M „). The gravitational field equations then assume 



once the Lorentz gauge condition W v \ v = is imposed. We are interested in 
the special retarded solution of (1) given by 



Setting T 00 = pc 2 and T° l = cj l , where p is the matter density and j is 
the matter current, inspection of solution (2) reveals that ft-oo = 4<!>/c 2 , h i = 



the form 



= - 




(1) 




(2) 
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—2Ai/c 2 , and hij = 0(c~ 4 ). Here $ and A have the interpretation of GEM 
scalar and vector potentials, respectively. Far from the source, e.g., <& ~ GM/r 
and A ~ GJ x x/cr 3 ; in general, <J> reduces to the Newtonian gravitational 
potential while A = 0(c _1 ). Neglecting terms of 0(c -4 ), we find that the 
Lorentz gauge condition reduces to 

19$ 1 

-- + V.(-A) = 0. (3) 
The spacetime metric is then of the form 

:^r)dt 2 - -(A • dx)cft + (1 + 2- 
tr c c 

The GEM fields are defined by 



ds 2 = -c 2 (l - 2-^)dt 2 - -(A • dx)dt + (1 + 2— )8 l0 dx l dx j . (4) 



E = -V$-~(^A) , B=Vx A, (5) 
which immediately imply the source-free GEM field equations 

VxE=-i|(iB),V.(iB) = 0, (6) 

while the gravitational field equations (1) imply 

„^ . „ , 1„ N 1 9„ 4irG . 

V-E = 47rGp, Vx(-B) = --E+— j. (7) 

These equations contain the conservation law for mass current dp/dt+V j = 0, 
as they should. Moreover, the gravitomagnetic field appears in (6) -(7) with a 
factor of ^ since qs/qE = 2. 

Let us recall that the gravitational potentials h^ u are gauge dependent; 
that is, under a coordinate transformation x^ — > a/^ 1 = rr^ — e^ 1 , — > ft^ = 
fyui/ + + to linear order. Therefore, 

and the Lorentz gauge condition is preserved provided De^ = 0. In general, the 
gauge freedom leaves the curvature tensor invariant but not the connection; 
however, it is possible to restrict the coordinate freedom in such a way that the 
GEM fields remain invariant. To this end, let e° = 0(c~ 3 ) and e 1 = 0(c~ 4 ); 
then, e", Q = 0(c~ 4 ) and with * = c 2 e°/4 we find that GEM potentials are 
defined up to the gauge transformation 

= $ - ~* , A' = A + 2V*, (9) 
c at 
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where * is a solution of = ("Lorentz gauge"). The GEM fields, which 
are in fact elements of the Christoffel connection, remain invariant under the 
gauge transformation in this approximation scheme. 

The motion of a test particle in the GEM field follows from the Lagrangian 
L = —m c ds/dt that can be evaluated to linear order in $ and A from (4). 
The result is 

V 1 V QtTTl 

L = -m.c 2 {\ - —)2 + 7717(1 + — )$ - — 7V • A, (10) 

where 7 is the Lorentz factor. If dA/dt = 0, we find that to lowest order in 
v/c, 

F = ?J! E| ?B -xB (11) 

c 

in analogy with the Lorentz force law. This follows from the fact that the 
canonical momentum of a test particle is given by p + (— 2m/c)A, where p is 
the kinetic momentum. The exchange of energy and momentum between the 
particles and the GEM field could naturally lead to the construction of the 
Maxwell stress-energy tensor for the GEM field. 

For the sake of convenience, it is possible to express the GEM potentials 
using the 4-vector notation as 

^" = (2*, A) (12) 

so that the Lorentz gauge (3) can be written as A^,^ = 0. We set c = 1 
in the rest of this paper, except where indicated otherwise, and define the 
GEM Faraday tensor as = A v ^ - A^ v . Then the GEM Maxwell 

equations are F^a] = an d , v — iirGj^, where = (p,j) is the mass 
current. It is then tempting to construct the analog of the Maxwell tensor as 
in electrodynamics, 

G T af3 = ^{T a 1 T p '< - \jf Si T 1& T^). (13) 

47T 4 

However, the physical significance of this quantity is doubtful as first pointed 
out by Maxwell [13] in his fundamental work on the dynamical theory of the 
electromagnetic field. The basis for Maxwell's considerations was the notion 
that the attractive nature of gravity would lead to a negative energy density 
for the field, while the electromagnetic analogy would imply a positive result. 
In fact, in our approximation scheme we can use instead the standard pseu- 
dotensor of general relativity that gives a negative energy density. This is 
the subject of section 4. 
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Finally, we digress here to point out the usefulness of the 4- vector notation 
introduced in (12) by considering the propagation of test electromagnetic radi- 
ation in the GEM field. In the absence of and in the background inertial 
spacetime, let us introduce the Kramers vector f = e + ib, where e and b are 
the electric and magnetic fields, respectively. Then Maxwell's equations for 
the radiation field can be expressed as V • f = and V x f = i df/dt. The 
latter equations can be written in the "Dirac" form using p a = —ihd a . If we 
now consider test electromagnetic radiation propagating in the GEM field by 
the standard substitution p a — > p a — qA ai where qe = —Tioj and qs — —2Tiuj 
for radiation of frequency u, then we recover the wave equation 

(V + 2iujA) x f = wNf, (14) 

where N — 1 + 2$ has the interpretation of the index of refraction [14]. It 
is interesting to note that in the Lagrangian (10) as v — > c the coefficient of 
$ becomes twice the energy of the particle, which is also the coefficient of the 
term v • A in (10); this is consistent with the relativistic substitution for the 
4-momentum p a — > p a — qA a with A a defined by (12). 

3 Gravitational Larmor Theorem 

In electrodynamics, the Larmor theorem has played an important role in the 
description of the motion of charged spinning particles. Explicitly, for slowly 
varying fields and to linear order in v/c and field strength, the electromagnetic 
field can be replaced by an accelerated system with translational acceleration 
&l = — qE/m and rotational frequency ujl = <?B/(2mc). For all charged par- 
ticles with the same q/m, the electromagnetic forces are locally the same as 
inertial forces; this is reminiscent of the principle of equivalence. However, 
q/m is not the same for all charged particles; hence, a geometric theory of 
electromagnetism analogous to general relativity is impossible. The geomet- 
ric treatment of gravitation is a direct consequence of the universality of the 
gravitational interaction. 

The gravitational Larmor theorem expresses Einstein's principle of equiv- 
alence within the GEM framework. To see this clearly let us consider an 
event in spacetime with GEM fields E and B. The corresponding Larmor 
quantities would be a L = —q E ~E/m = E and uj l = q B H/(2mc) = — B/c, in 
accordance with our convention. In the neighborhood of this event, we can 
replace the GEM field by a neighborhood in Minkowski spacetime that has 
translational acceleration a^ — E and rotational frequency lol = B ( with 
c = 1). The geodesic coordinates in the neighborhood of such an accelerated 
system would be X' 1 = (r, X) such that the Minkowski metric can be written 
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as ds 2 = g llv dX^dX v with 

5oo = -(1 + a L • X) 2 + (oj l x X) 2 , 



(15) 



9oi = (wl x X)i , gjj = 



(16) 



Let us note that to linear order in X, the metric (15)-(16) would become the 
same as (4) once we set $ = — az, • X and A = x X as in section 2 

and neglect spatial curvature. Then the GEM fields corresponding to these 
potentials are E = — V<& = and B = V x A = -w^, thus confirming the 
validity of the gravitational Larmor theorem in our approach. 

Einstein's heuristic principle of equivalence is traditionally stated in terms 
of the translational acceleration of the "Einstein elevator"; however, it is clear 
from the analysis of this section that a Larmor rotation of the elevator is nec- 
essary as well to take due account of the gravitomagnetic field. 



4 Stress-Energy Pseudotensor 

The standard Landau-Lifshitz pseudotensor is useful in a global sense for 
asymptotically flat spacetimes. It is therefore appropriate in our linear per- 
turbation analysis around Minkowski spacetime for r » 2GM/c 2 ; however, we 
use t^y to compute the local stress-energy content of the GEM fields. We find 
that 



16ttG t' 



-UE 2 + Y,(Aij + A jti )(A id + Aj.i) - 6($, f 



-UE*A h0 - -A\ A ifi , 



(17) 



4ttG t 0i 




(18) 



WirG tii 



- U l3 E 2 ) + 4(B i S j + h tJ B 2 ) 
-6[EiA jfi + E 3 A lfi + <%($,o f - 5 i3 E k A kfi ] 
-7{A ifi A jfi + ^6ijA k fi A kfi ). 



(19) 
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For a stationary source with angular momentum J = Jz, the dominant 
contributions are 



nn 7 GM 2 , n ,. 1 GMJ . n 7 
87r r 4 Z7r r b 

1 CM 2 

t ij = — H_a;V. (21) 

It is important to note that the field energy density is negative and that there 
is a flux of energy that circulates around the mass in a direction opposite to 
its sense of rotation. The flow velocity, v l g = t 0l /t 00 , is given by 

Vs = lW s&n9 ^ (22) 

For the sake of comparison, it is interesting to point out that the Maxwell 
tensor (13) would give the same result as (22) except that the numerical coef- 
ficient would be 1 instead of |. 



5 GEM Stress-Energy Tensor 

Thus far our perturbative treatment of GEM has involved a fixed global back- 
ground inertial frame. This is contrary to the spirit of the general theory of 
relativity; therefore, it is important to develop an alternative approach for an 
arbitrary gravitational field. To this end, imagine a congruence of test geodesic 
observers in a gravitational field. Choosing an observer in the congruence as 
a reference observer, one can set up a Fermi coordinate system in its neigh- 
borhood. Instead of a global background inertial frame, we thus have a local 
inertial frame along the worldline of the fiducial observer. The spacetime met- 
ric in the Fermi system is given by 

goo = -l-R ioj(T)X i X* + ..., (23) 
g 0l = -^Ro jik (T)Xi x k + (24) 

9ij = S t j - -R lk jiX k X l + (25) 
where X^ = (r, X) are Fermi coordinates and 

Raf3~/5 = -^wpctA^A^A^jA^) (26) 
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are the components of the Riemann tensor projected on the Fermi-propagated 
orthonormal tetrad frame s (r) of the fiducial observer at the origin of spatial 
Fermi coordinates (r, 0). Let us note that the deviations from the local fiducial 
inertial frame arc of the form 

$(t,X) = ~R <HOj (T)X i X' + ..., (27) 

Mt,X) = l -R 0jlk {T)Xi X k + (28) 

where in this GEM approach we simply ignore the spatial curvature. The 
corresponding GEM fields would then be 

£?i(r,X) = R 0i0j (T)Xi + (29) 

Bi(r,X) = ~e ijk R jk0 i{T)X l + (30) 

where we have used our previous conventions. 

From the viewpoint of our reference observer, the other test observers in 
the congruence move in accordance with the generalized Jacobi equation [15] 



^4- + R 0i0j X j + 2R ikj0 V k X j + {2R 0kj0 V i V k 

+^R ikj iV k V l + ^R ok] iV l V k V l )X^ = 0j (31) 

which is valid to first order in the separation X while the rate of separation 
V = dX-fdr could be arbitrary (|V| < 1). Restricting the deviation equation 
(31) to first order in the relative velocity as well, we find that it can be written 
in the form 

d 2 X. 

m— = q E E + q B V X B, (32) 

with qs = —m and qs = —2m as before. It is interesting to note that the 
gravitoelectric and gravitomagnetic fields (29)- (30) are given by the electric 
and magnetic components of the Riemann tensor, respectively. 

It is possible to combine (29) and (30) and define a GEM Faraday tensor by 
F a p = —R a poiX l to linear order in X. Here F$i = —Ei and = £ij k B k as in 
standard electrodynamics; therefore, this treatment is in general different from 
our perturbative approach. This Faraday tensor satisfies Maxwell's equations 
F[a/3,j] = and F a ^,p— AirJ a with 47rJ Q ,(r, 0) = — R 0a along the fiducial 
worldline. Therefore, the classical field theory of the GEM field in the Fermi 
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frame can be developed in the standard manner [16]. We are particularly 
interested in the Maxwell stress-energy tensor T a/3 associated with F a p-. this 
is given, in a similar form as (13), by 

G T^(r, X) = -L (i r 70j i?%. - -^R^R^^X 1 ^. (33) 

The fact that T a/3 vanishes for the reference observer is consistent with Ein- 
stein's principle of equivalence. Away from this fiducial observer, T a @ is 
nonzero in general; but, we could have just as well chosen a different ref- 
erence observer in the congruence and then the situation would have been 
reversed. In addition, we recall that physical measurements of field energy and 
momentum cannot be performed at a point and require an averaging process 
as already emphasized by Bohr and Roscnfeld [17]. Therefore, we average (33) 
at an event (r, 0) along the reference trajectory over an infinitesimal sphere 
of radius e£. Here < e <C 1 and I is an invariant length scale (e.g., i could 
be GM/c 2 or the Planck length when a natural length scale is absent in the 
case under consideration). Then < X % X^ >= Co(ei) 2 Sij, where Co = | or 
■| for averaging over the volume or the surface of the sphere, respectively; in 
any case, this constant coefficient can always be absorbed in the definition of 
I. Hence 

< T afj >= 47rG T^pffA^Aj^A^A^), (34) 
where T puprT (x) is the Bel tensor given by 

Tfj, upcT (x) = -(R^pqR^ ^ + R^e^R^ p C) - -^g liV R a { ipi R a(i CT 7 . (35) 

This tensor was first defined by Bel in 1958 for Einstein spaces on the basis of 
a certain analogy with the Maxwell stress-energy tensor [18]. It is symmetric 
in its first and second pair of indices and traceless in the first pair. We define 
the GEM stress-energy tensor to be T^x), where the stress-energy measured 
by an arbitrary observer is essentially given by the projection of the Bel tensor 
on the tetrad of the observer, i.e. 

P ~ » v p 

T (")(/3) = Q T ^P<? X (a) X (P) X (0) X (0) ■ ( 36 ) 

For a Ricci-flat spacetime, the curvature tensor reduces to the Weyl con- 
formal tensor C^ pa and the Bel tensor reduces to the totally symmetric and 
traceless Bel-Robinson tensor T^ prT given by 

T^upa = — {C ^p^Cv^ ^ + C^ a ^C v ^ p^) — — gpvgpaC a p~ f sC a P~ <s . (37) 
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In this case, the spatial components of the curvature tensor (that were ignored 
in our construction of the GEM stress-energy tensor) are in fact given (up to 
a sign) by its electric components. Therefore, T^ u (x) reduces in this case to 
the gravitational stress-energy tensor T* ^ v given by 

e v 

T* (a)(0) = q T ^P<T^( a )^(0)^(O)^(O)- ( 38 ) 

This GEM derivation of the physical content of the Bel and Bel- Robinson 
tensors has made it possible to provide new definitions for the average gravita- 
tional stresses (such as the pressure of gravitational radiation) up to a positive 
multiplicative factor [16]. 



6 Gravitational Poynting Flux 

The GEM approach based on the curvature tensor can be used in arbitrary 
gravitational fields. A simple application of the gravitational stress-energy 
tensor to the exterior field of a rotating source implies the existence of a steady 
flow of gravitational energy circling the source. This situation is analogous 
to the electromagnetic case of a charged rotating mass; that is, there is a 
steady Poynting flux of electromagnetic energy circulating around the mass. 
In connection with the physical interpretation of the stress-energy tensor, Pauli 
has noted that in a purely electrostatic (or, alternatively, magnetostatic) field, 
there is no momentum density, but there is momentum current [19], while 
Feynman has discussed some counter-intuitive features of the Poynting flux 
[20]. Nevertheless, one can study the general physical characteristics of this 
gravitational energy flow using the GEM stress-energy tensor [16]. 

For the exterior of a rotating mass (r 3> 2GM/c 2 ), the dominant contri- 
butions to the Bel-Robinson tensor are 

T( 0)(0) = ^ )T (0) W = 9^ ax . )i) (39) 

T (i)U) = ^Ml {2Sij -3x i ^), (40) 

for freely falling observers. 

The three approaches to the stress-energy content of GEM have this feature 
in common: They all predict that the gravitational energy will steadily flow 
around a rotating mass with a velocity 
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where A is a numerical factor equal to 1, y, or 3 for the three cases we have 
discussed. It is therefore interesting to examine certain general features of this 
flow velocity. We find that V-v g = 0, so the flow is divergence-free; in fact, the 
streamlines are circles in planes perpendicular to J. Moreover, the vorticity of 
the flow, ojg = V x v g = 2A( J /Mr 2 ) cos 9 f , vanishes in the equatorial plane 
and is maximum along the rotation axis. Let us note that v 9 = Vx ty g , where 
tyg is the gravitational stream function given by ^ g = X(J/M) cos 8 f . In fact 
Lo g = — V 2 \& g , since V • ^ g = 0. The circulation of v s , 

C g = jw g -d\, (42) 

around a streamline is given by 2ir\(J/M) sm 2 6; therefore, the circulation is 
independent of the radial distance r and vanishes on the axis of rotation. In 
the equatorial plane, the circulation C' g is simply proportional to the specific 
angular momentum of the source, just like the gravitomagnetic clock effect. It 
would be interesting to explore the observational aspects of the steady gravi- 
tational energy flux around rotating masses. 
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